Quantum pumping and dissipation: from closed to open systems 
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Current can be pumped through a closed system by changing parameters (or fields) in time. The 
Kubo formula allows to distinguish between dissipative and non-dissipative contributions to the 
current. We obtain a Green function expression and an S matrix formula for the associated terms in 
the generalized conductance matrix: the "geometric magnetism" term that corresponds to adiabatic 
transport; and the "Fermi golden rule" term which is responsible to the irreversible absorption of 
energy. We explain the subtle limit of an infinite system, and demonstrate the consistency with the 
formulas by Landauer and Buttiker, Pretre and Thomas. We also discuss the generalization of the 
fiuctuation-dissipation relation, and the implications of the Onsager reciprocity. 



Linear response theory (LRT) is the traditional theo- 
retical tool for dealing with the response of driven sys- 
tems [H H, H, Q . It offers an expression (the Kubo for- 
mula) for the generalized susceptibility, and hence for 
the generalized conductance matrix. It has been realized 
that in the adiabatic limit the Kubo formula reduces to 
an expression for "geometric magnetism" j^. In case of 
electrical current calculation the latter gives the "adia- 
batic transport" of charge Outside of the adiabatic 
regime the response includes an additional "dissipation" 
term 8]. The latter determines the rate of irreversible 
energy absorption, which is caused by Fermi-golden-rule 
transitions between energy levels. 

Recently there was much interest in analyzing the re- 
sponse of open systems that are connected to reservoirs. 
The analysis has been based on the S'-matrix formalism, 
leading to the Landauer formula 0, 13 , and more gener- 
ally to the Biittiker, Pretre and Thomas (BPT) formula 
0. A major motivation for the present work is the real- 
ization that the relation between the BPT formula and 
the Kubo formula has not been clarified. In particular 
the notion of "adiabatic pumping" in the context of an 
open system has been left obscured, and some confusion 
has arose regarding the role of dissipation in the pumping 
process HEElnl. 

The purpose of the present work is to analyze the re- 
sponse of closed isolated systems [Tllll2j|. and in particu- 
lar to consider the special limit of an infinite system (no 
reservoirs!). Thus we are going to construct a bridge be- 
tween the LRT formulation and the BPT formula. This 
is of great practical importance, because the assumed 
open geometry of the S'-matrix formulation is in many 
cases an idealization. It is clear that Kubo formula al- 
lows a straightforward incorporation of finite-size, exter- 
nal noise, environmental and possibly also many-body ef- 
fects. A major step in constructing this bridge, had been 
taken in Ref.^^i where the authors start with the Kubo 
formula for the electrical conductivity and end up with 
the Landauer formula which relates the conductance to 
the transmission of the device. We are going to see that 
the general case, which deals with the generalized con- 
ductance matrix, and hence incorporates adiabatic trans- 
port, is much more subtle. 

Consider a closed isolated system whose Hamiltonian TL 



depends on several control parameters Xj . An example is 
presented in Fig. 1, where xi and X2 are gate voltages and 
X3 is the magnetic flux through the loop. The general- 
ized forces are conventionally defined as — —dTC/dxk- 
Note that is the electrical current. In LRT the 
first order contribution to {F'') is related to Xj{t) by 
a causal response kernel a''^ {t ~t'). The Kubo expres- 
sion for this response kernel is a^^ {t) — Q{t) K^^{t), 
where K''i{T) = {i/TL){[F''{T),F^ {0)]), and Q{t) is the 
step function. The Fourier transform of a^^ (r) is the 
generalized susceptibility ^''^ (w). The generalized con- 
ductance matrix is: 



Qkj 



lim 



Im[x'^^(w)] 



K^^{T)TdT (1) 



Thus in the limit of zero frequency the non-trivial part 
of the response can be written as a generalized Ohm law 



(2) 



where following 5] the generalized conductance matrix is 
written as a sum of a symmetric matrix r] 



kj 



V 



,ik 



that 



represents the dissipative response, and an antisymmetric 
matrix B^^ — —B^^ that represents the non-dissipative 
response (also called "geometric magnetism"). 

For a device as in Fig.l, and zero temperature occupa- 
tion of non-interacting (spinless) Fermions, we find below 
that the dissipative part of the response is 



- trace Im[G+] F^ Im[G+] 
dS^ dS ' 
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dxi dxj 
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where G = 1/{E—Ti±i0) are Green functions of the cor- 



responding open system, and Im[G 



+1 - 



-zi(G+-G-). 



For the non-dissipative part of the response we find 
ih 
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where the second equahty holds for fc = 3, and allows the 
determinatfon of the electrical current (F^) via a speci- 
fied lead A. The last term is defined in Eq. H26(l . The pro- 
jector Pa restricts the trace operation to be over the spec- 
ified lead channels. In the absence of magnetic field the 
remaining component is B^^ — 0, while r]^^ — rf^ = 
as expected from the Onsager reciprocity relations (see 
last paragraph). Disregarding the last term in Eq.®, 
the sum of © and l@J for fc = 3 coincides with the BPT 
formula, which can be written in our notations as: 



G3.-^^trace fp.|^5t 



fBPTl 



(7) 



We show later that this reduces for j = 3 to the Landauer 
formula which relates the electrical conductance G^^ to 
the transmission of the device. 

Below we explain how to derive the expressions for r}''^ 
and B'^^ starting from the Kubo formula Eq.(^). Later 
we discuss further physical implications of our results. 
Assuming zero temperature Fermi occupation up to en- 
ergy Ep , standard textbook procedure 0, 0, |M 01 leads 
to 



^f^J2^nm KEf - Era) P,i„ S^Ep - E,,) (8) 



where the overline indicates that the delta functions are 
smeared. If the system were not isolated, the "broad- 
ening" r of the energy levels would be determined by 
the interaction with the external environment j^]. But 
we assume a closed isolated system. Still we argue Q 
that in case of a quantized chaotic system the levels ac- 
quire an effective width F — ((?icrF/A^)|a;|)^/'^A, where 
A is the mean level spacing, and ap is the root-mean- 
square value of the near-diagonal matrix elements (see 
remark |14|). Therefore we find two possibilities: In the 
adiabatic regime (F <C A) the dissipative conductance is 
zero (rj = 0), while in the non- adiabatic regime (F > A) 
the dissipative conductance acquires a well defined finite 
value, which is not sensitive to F, and can be calculated 
using Eq.|(31l. A similar claim holds regarding B*^-' , but 
the details are much more subtle: We start with the stan- 
dard expression ^| 



r=o 



Ofc --/p ^ \ " Im [-F^m-P'mn] 



(9) 



where f{E) is the Fermi occupation function (later we 
take the limit of zero temperature). Incorporating F, 
and exploiting the antisymmetry of the numerator with 
respect to n <^ m interchange we get 



E 



nm mn 

(p„-p„)'+(r/2)' 



;{fiEnyf{Ern)) (10) 



The numerator, on the average, depends mainly on the 
difference r — m — n, and it is non-negligible within a 
bandwidth \Em ~ En\ < A^. We further discuss the 



bandwidth issue in the next paragraph, and explain that 
in the limit of a very long wire A ^ F ^ A(,. This 
means that in this limit F serves like the infinitesimal 
iO in the definition of the Green fimctions G*. Con- 
sequently, the sum in Ea. (|10l) . which is of the form 
T,n,,n9{n~m){f{En)-f{E,n)) = Er^gir), leads after 
some straightforward algebra to Eq.Q. 

Formally there is an optional derivation that leads to 
(O and (0. The kernel K^^{t) is related to the sym- 
metrized correlation function C"-'(r) = {^{F^{t)F^ (0) + 
F^ (O)F' (t)) . The quantum mechanical derivation of this 
subtle relation is discussed in Appendix D of 0. If we 
use this relation we get from Eq. 1^ an extremely simple 
(and useful) result: 



(11) 



which can be regarded as the generalization of the fluctu- 
ation dissipation relation. The fluctuations are described 
by C'^^ (uj) which is defined as the Fourier transform of 
C*^ (r) . It follows from this definition that 



2Trh 



C'^^u;) = — FH^FLr. 



(12) 



For the device of Fig. 1 the mean level spacing is A cx l/L 
where L is the length of the wire. The above relation im- 
plies that the bandwidth of the mn matrix is Af, ~ h/rd, 
where the classical correlation time Td is determined by 
the chaotic motion inside the dot. It is also clear that 
C^^{lo) oc l/L, and therefore a\ cx (1/L)^. Hence we 
get that F cx (l/L)^/-^, implying that the limit L — > oo 
(keeping constant Fermi energy) is non-adiabatic, and 
that A ^ F ^ Af). Assuming for simplicity that there is 
no magnetic field, one easily derives the expressions 



1 

A 



UJ 



duj 
2^ 



(13) 

for j==l,2 (14) 



which are equivalent to those obtained in the previous 
paragraph. Note that C^^{t) with j — 1,2 is antisym- 
metric with respect to r, and therefore —iC^^{uj)lbJ is a 
real symmetric function. 

We turn back to the formal derivation. We want to get 
exact expressions for the elements of the conductance ma- 
trix, for the device of Fig. 1, in the non-adiabatic limit 
of large L. The location of the particle is specified by 
r ~ (r, s), where r is the coordinate along the ring, and 
s is a transverse coordinate. Optionally we can spec- 
ify the location along a lead using a radial coordinate 
r, while the surface coordinate s distinguishes different 
points that have the same r. We shall refer to r = as 
the boundary of the scattering region. The channel basis 
is defined as (r, s\a, r) = Xa{s) 5(r— ra(r)), where a is the 
channel index. The wavefunction in the lead regions can 
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be expanded as follows: 



ikar\ 



(15) 



Following we define an operator which can be iden- 
tified with the imaginary part of the self energy of the 
interaction of the dot with the leads: 

f ^ |a, 0)hva{a, 0| = S{r) (g, ^ \a)hva{a\ (16) 



where Va = (?ifca/mass) is the velocity in channel a. The 
matrix elements of the second term in Ea. (|16|l are 



f (S, s') = ^ Xa{s) flVa Xa(s') 



(17) 



Using standard procedure (see section 3.4 in the 
Green function in the leads, inside the scattering region 
(r < 0), can be expressed using the S matrix: 

G+(r,5|0,s')= (18) 

where k = diagjfca} is a diagonal matrix. Now we are 
fully equipped to convert Eq. |j3|into an S'- matrix expres- 
sion. Using the identities (see Q for Ea. H19|l 'l 



Im[G+] = -iG+EG- = -iG-FG^ 



dxi 



we obtain 



— trace 
An 



dG+^dG-^ 



(19) 
(20) 

(21) 



Using the definition of f and Eq. ifT^ we get Eq. Q . 

The derivation of the S matrix expression Eq.lj^Jl for 
B''^ is much more subtle, and requires a preliminary dis- 
cussion of the definition of the current operator. Consider 
a ring geometry, and assume that the current is driven 
by the flux In order to have a better defined model we 
should specify what is the vector potential A{r) along the 
ring. We can regard the values of A at different points in 
space as independent parameters (think of tight binding 
model). Their sum (meaning f A{r)-dr) should be $. So 
we have to know how $ is "distributed" along the ring. 
This is not just a matter of "gauge choice" because the 
electric field £{r) = —A{r) is a gauge invariant quantity. 
The transformation A i-^ A + VA(r) for a time depen- 
dent filed is not merely a gauge change: A gauge transfor- 
mation of time dependent field requires a compensating 
replacement of the scalar potential, which is not the case 
here. So let us define a flux <I>^ which is associated with 
a vector potential that is concentrated across a section 
r = r^i of a given lead. For the later derivation it is es- 
sential to assume that the section r = is contained 



within the scattering region (see Fig. 1). The generalized 
force which is associated with is — Ia, the current 
through this section. Namely 



Ia = - 



dn 



= ^e{v S{r ^ ta) + d{r - rA)v) (22) 
= ie/h)[tAP+ -TaP-] (23) 



where v is the r component of the velocity operator. The 
last equality involves new definitions. We define 



aeA 



\a,rA)hva{a,rA\ 



(24) 



We also define projectors P+ and that project out 
of the lead wavefunction Ea. p5|l the outgoing and the 
ingoing parts respectively. These projectors commute 
with Fa- Furthermore, note that P~^G~^ — G"*", and that 
P^G+ = 0, and that G^P^ = etc. Using these extra 
identities one obtains the following expression: 



A-Kl 



trace 



f.^fG^ 

OX-i 



(25) 



Using the definitions of f and Fa, together with Ea. (|18|l . 
followed by a straightforward algebraic manipulation, one 
arrives at Eq.® with the additional term 



= —Re 



dS 
' 9a;,- 



trace I Pa ^e'^'^''* 



(26) 



where Pa is a projector that restrict the trace operation 
to the aeA lead channels. 

For the simple ring geometry of Fig. 1, we have a left 
lead [h e B) and a right lead (a g A) channels, and the 
S'-matrix can be written as 



S = 



Pa = 




1 



where <j) — e^j^/h. Using the identity 



= iUPaS-SP,) 
n 



(27) 



(28) 



one can derive the relation that has been stated between 
Eq.Q, Eq.® and the BPT formula Eq.Q. Further- 
more, assuming that there is an electro-motive force — $b 
which is induced in the other lead, one obtains from BPT 
^33 _ (e^/27r?i)trace(tABiAB) which is the Landauer for- 
mula. The application of this procedure to multi-lead 
systems is a straightforward generalization. 

For an open system it is evident that the current which 
is emitted (say) through the right lead, does not have to 
be equal to the current which is absorbed by the left 
lead. The reason is that charge can be accumulated in 
the dot region. But for a pumping cycle one realizes 
that the integrated current (pumped charge) is a well 
defined (lead independent) quantity. Similar observation 
holds in case of a closed system. Assume for example 
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that the left lead is blocked. In such case raising the dot 
potential will cause an emission of charge through the 
right lead, while the current through the left lead is zero. 
The emitted charge is accumulated in the "wire". But 
for a full cycle the original charge distribution is restored, 
and therefore the integrated charge (Q) becomes a well 
defined (lead independent) quantity. The additional term 
Ea. (|26|l gives a zero net contribution for a full pumping 
cycle. This term implies that the current is not uniform 
within the lead. The current has a modulation in the 
radial direction (r), with a spatial period that equals half 
the De-Broglie wavelength at the Fermi energy. This 
reflects that the net transported current corresponds to 
translation of a standing wave which is associated with 
the last occupied level. 

More subtle is the value of Q for a full driving cycle. 
In contrast to a previous wrong statement we have 
argued pTl| that for a strictly adiabatic driving cycle, in 
the absence of magnetic field, the transported charge Q 
is at best approximately quantized (say Q w 1 in units of 
the elementary charge). The deviation is related to the 
Thouless conductance of the device, and can be either 
positive or negative In contrast to that, with the 

BPT formula the correction to Q « 1 is always nega- 
tive. On the basis of our derivation we can conclude the 
following: The deviation from quantization in a strictly 
adiabatic cycle is related to the contribution of the neigh- 
boring level. If the degeneracy with this level is located 
in the plane {xi,X2tX^=Q) of the encircling cycle, then 
the correction is positive. If the encircled degeneracy is 
off plane, then the correction is negative. The effect of 
non-adiabaticity (F > A) is to screen the contribution 



of the neighboring levels, which is the reason for having 
always a negative correction from the BPT formula. 

The role that dissipation may have in pumping is re- 
stricted, merely by the realization that the BPT for- 
mula is related to the Kubo formula. The Onsager reci- 
procity relation imply that in the absence of magnetic 
field the conductance matrix G^^ should be symmetric 
(antisymmetric) with respect to the permutation of the 
indexes {k,j), depending on whether F'' and trans- 
form (not) in the same way under time reversal. This 
means that shape deformations lead to dissipation via ry*-' 
with ij <3, while the electrical current is determined ex- 
clusively by the non-dissipative terms B^^ and B^^. This 
should be contrasted with the response to electro-motive 
force which is purely dissipative: Both the current and 
the dissipation are exclusively determined by the Ohmic 
conductance rj'^^. Thus, in the absence of magnetic field, 
we have a clear cut distinction between the dissipative 
and the non-dissipative contributions to the response. 

In summary, starting with the Kubo formalism, we 
were able to find expressions for the dissipative and for 
the non-dissipative parts of the response, and to illumi- 
nate the role of non-adiabaticity in the limiting case of an 
infinite system. In contradiction with past speculations, 
we were able to demonstrate that the switch to an open 
system does not necessitate an extra dissipative term. 
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FIGl. Illustration of a closed system. The dot potential is con- 
trolled by gate voltages xi and X2- The flux through the loop is 
a;3=<t>. The scattering region {r<0) is represented by an S matrix. 
The length (L) of the wire is assumed to be very large. 
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NOT PART OF THE PAPER 



Going from Eq.(lO) to Eq.(5) 

In this appendix we give the "straightforward algebra" that leads from Eq.(lO) to Eq.(5). 

Within the framework of the approximation which is discussed in the text this sum is of the form 

B^'i = Y,9{n-m)x{f{En)-f{Em)) 

Using the notation r = m — n it leads to 

r 

We make the substitutions 

E — E 



9{r) = 



A 



nm mn 



where n is the index of an arbitrary energy level in the vicinity of the Fermi energy {En ^ Ep)- 
Now we have the expression 



B^i = -in— V E^-En J 



Thanks to V we can make the replacement E^ i— > Ep. 
Next we can use (in reverse) the approximation 



Y,{n\A\n)5{En-EF) « \{n\A\n) 



where the overline indicates a smeared delta function. 
Hence we get 

As explained in the text, in the limit of a very long wire F is like the infinitesimal iQ. 
Hence we get Eq.(5). 
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APPENDIX D OF REF.[8] 



Expressing K{ij) using C{uj) 

We can use the following manipulation in order to relate K'^^{lo) to C*^ (tj), 

k^^iiu) = ^/(i?„) (29) 

n 

nm 



nm 



nm 



5] /'(£;„) H 



where we use the notation ujnm = {En — Em)/Ji- The third line differs from the second line by permutation of the 
dummy summation indexes, while the fourth line is the sum of the second and the third lines divided by 2. In the last 
equality we assume small u. If the levels are very dense, then we can replace the summation by integration, leading 
to the relation: 

k'Hu;) = I g{E)dE f{E) k%{u,) = ~iu ( g{E)dE f'{E) C%{u,) (30) 



where k^{uj) and C^(w) are microcanonically smoothed functions. Since this equality hold for any smooth f{E), it 
follows that the following relation holds (in the limit ^ 0): 



^l(-) = (31) 

If we do not assume small oj, but instead assume canonical state, then a variation on the last steps in Ea. l29|l . using 
the fact that {f (En)— f (Em)) / if {En)+ f{Em)) — tanh((i5„— £^m)/(2T)) is an odd function, leads to the relation 

Kffi{cu) = X -5- tanh f — ) Cff! (cu) (32) 



Upon substitution of the above expressions in the Kubo formula for G'-' , one obtains the Fluctuation-Dissipation 
relation. 



Note added: For a low temperature Fermi occupation Ea. H3U|) can be written as 

k^{u) = zu: giEp) C%^{u;) (33) 

Formally this is valid only if hu) ^ T. But in fact if T <C Ep the result should not be sensitive to Ep- Therefore it 
can be argued that this relation holds globally, and we can set T — Q. It follows that 

K''{r) = ^^g[EF)C%^{T) (34) 

This can be used in order to derive Eq.(ll) for Eq.(l). We see that it is a subtle relation which should not be regarded 
as a trivial identity. 



